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Abstract.  Vibrational  and  electron  kinetics,  as  well  as  electronic  excited  states  and  ionization  kinetics,  have  been 
coupled  self  consistently  solving  master  and  Boltzmann  equation  at  the  same  time.  These  kinetics  have  been 
inserted  in  a  quasi  one-dimensional  code  for  the  high  enthalpy  nozzle  expansion. 


INTRODUCTION 

The  supersonic  expansion  through  a  nozzle  flow  is  characterised  by  strong  non-equilibrium 
conditions.  Multitemperature  approach  does  not  describe  correctly  the  vibrational  distributions.  In 
previous  works  [1-3]  we  have  shown  the  non-Boltzmann  character  of  the  vibrational  distributions  in 
N2/N  [1],  and  air  mixtures  [2,3],  especially  when  atom  recombination  is  dominant;  starting  with  very 
high  temperature  in  the  reservoir,  the  gas  is  partially  dissociated  and  during  the  expansion  the 
recombination  is  strong,  generating  non-Boltzmann  vibrational  distributions.  As  a  consequence,  the 
global  dissociation  rates  do  not  follow  the  Arrhenius  dependence  on  the  local  temperature. 

If  the  reservoir  temperature  is  high  enough  to  have  ionization  degree  of  the  order  of  lO'^-hlO'^, 
electron-molecule  (e-M)  collisions  cannot  be  neglected.  The  conditions  of  the  system  during  the 
supersonic  expansion  can  be  compared  to  post-discharge  [4]:  energy  trapped  in  internal  degrees  of 
freedom  (vibrational  and  electronic  excited  states)  and  electron  kinetic  energy  are  very  high,  therefore 
the  system  should  relax  to  reach  the  equilibrium  with  the  translational  temperature  of  molecules.  In  the 
nozzle  supersonic  expansions  there  is  a  further  difficulty  due  to  the  fact  that  gas  temperature  and 
pressure  change  during  the  expansion,  while  these  quantities  are  constant  in  post-discharge  systems. 

Even  if  ionization  degree  is  high,  electron-electron  collisions,  though  very  important,  are  not 
sufficient  to  overcome  the  effect  of  superelastic  and  inelastic  collisions,  and  therefore  non-Maxwell 
electron  energy  distributions  are  obtained  [4].  To  model  non-equilibrium  vibrational  distributions  in 
presence  of  non  maxwellian  electron  energy  distributions  it  is  necessary  to  solve  contemporary  the 
master  equations  for  the  vibrational  kinetics  and  the  Boltzmann  equation  for  free  electrons  (self- 
consistent  model). 

The  self-consistent  model  has  been  applied  in  flow  conditions  for  the  boundary  layers  of  hypersonic 
re-entering  body  [5]  and  for  nozzle  supersonic  expansion  [6],  neglecting  ionization  and  electronically 
excited  state  kinetics,  showing  strong  effect  of  the  self-consistent  coupling  in  both  vibrational  and 
electron  energy  distribution  functions  (eedf);  electronically  excited  states  fraction  as  well  as  ionization 
degree  have  been  considered  constant.  In  recent  work  [7]  we  have  applied  the  model  to  the  conic  nozzle 
adding  the  electronically  excited  states  and  ionization/recombination  kinetics.  The  results  have  shown 
that  as  supposed  in  previous  works  ionization  degree  is  practically  constant,  because  the  pressure  is 
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very  low,  but  electronically  excited  states  change  in  a  range  of  many  orders  of  magnitude.  Strong  non- 
non-equilibrium  distributions  are  obtained,  with  a  strong  effect  of  superelastic  collisions. 

In  this  work  we  have  applied  the  model  used  in  [7]  to  F4  nozzle  [8]  to  study  the  effect  of  a  different 
nozzle  profile. 


NUMERICAL  MODEL 

The  model  used  to  calculate  the  nozzle  fluid  dynamics  is  based  on  the  quasi  one-dimensional  steady 
Euler  equations  [1-3] 
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for  reacting  flows,  where  "p^v"  is  the  mass  density  of  the  species  s  in  level  v.  These  equations  are 
completed  by  the  definition  of  the  enthalpy  for  a  non-equilibrium  flows 

h  =  hj  +  hj  +  hg  (2) 


as  the  sum  of  the  translational  (eq.  3),  internal  (eq.  4)  and  free  electron  (eq.  5)  enthalpy 


h^  =  RT 


aRT 


hi  = 


EXs(H{  +  f,^e^^ 
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(3) 

(4) 

(5) 


The  quantity  "f^y"  is  the  energy  level  distribution  or  the  fraction  of  heavy  particle  of  s-th  species  in  the 
v-th  level  and  "n(e)"  is  the  eedf. 

The  contribution  of  reactions  is  included  in  the  mass  continuity  equation  for  any  single  species  (see 
eq.  1)  where  the  source  term  can  be  calculated  from  the  zero-dimensional  kinetic  equations  obtaining  the 
following  expression  for  species  mass  fraction 
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This  model  has  been  applied  to  study  vibrational  and  chemical  kinetics  and  can  be  easily  extended  to 
the  electronically  excited  states.  Some  rearrangements  are  needed  to  extend  the  nozzle  model  when 
electron  kinetics  is  coupled  self-consistently  with  vibrational  and  excited  level  (vibrational  and 
electronic)  kinetics  (see  fig.  1). 


FIGURE  1.  Scheme  of  the  self-consistent  coupling  between,  heavy  particle,  electron  kinetics  and  Euler  equations. 

As  starting  point  we  consider  the  Boltzmann  equation  for  free  electrons  in  the  homogeneous  and 
almost  isotropic  (two  terms  in  spherical  harmonics  expansion)  approximation  [9].  The  time  dependent 
eedf  can  be  calculated  solving  the  equation 

dnfe)  a  ,  . 

^  =  +  (7) 

where  Jg,  Jgj,  Jgg  are  the  flux  along  the  electron  energy  axis  due  respectively  to  electric  field  (E), 
electron-heavy  particle  (el)  and  electron-electron  elastic  collisions  and  and  S5,„p  are  the  source  terms 
due  to  inelastic  (eq.  8)  and  superelastic  (eq,  9)  collisions  (i<j). 

e(£)  -I-  N2(i)  —>  e(£-eij)  -i-  N2(j)  (8) 

e(£)  -I-  N2(j)  — >  e(£-i-£jj)  +  N2(i)  (9) 

The  "J"  terms  depend  on  the  derivatives  (first  and  second)  of  the  eedf  while  "S"  terms  make  the 
electron  to  jump  from  one  energy  to  another  and  are  proportional  to  the  eedf.  All  the  terms  are  linear 
but  "Jgg"  which  is  quadratic  in  the  eedf  [11]. 

The  numerical  approach  to  the  solution  of  the  nozzle  equations  has  been  widely  described  in 
previous  works  [1-3,7].  For  N2  gas,  the  species  considered  are  N2,  N2‘'‘,  N,  N'*',  e".  Only  for  nitrogen 
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molecules,  internal  level  kinetics  has  been  taken  into  account,  considering  45  vibrational  levels  and  four 
electronically  excited  states  (N2(A)),  B^Ilg  (N2(B)),  C^Hg  (N2(C))  and  a'^E„'  (N2(a)). 

We  have  used  the  same  vibrational  and  dissociation/recombination  kinetics  used  in  [1-3] 


N2(v)-I-N2(w)  N2(v-1)-I-N2(w-I-1) 

(10) 

N2(v)+N2  N2(v-1)-I-N2 

(11) 

N2(v)-i-N  N2(w)-tN 

(12) 

N2(v)+N2(v')  ^  N2(v-1)-i-2N 

(13) 

N2(v')+N2  ^  2N-tN2 

(14) 

N2(v)-i-N  3N 

(15) 

To  this  model  we  have  added  electronically  excited  states  and  ionization/recombination  kinetics  [7] 


N2(A)-i-N2(A)  N2(B)-i-N2(v=8) 

(16) 

N2(A)-i-N2(A)  ^  N2(C)-hN2(v=2) 

(17) 

N2(A)+N2(v>6)  N2(B)-i-N2(v-6) 

(18) 

N2(A)+N  N2(v<10)-i-N 

(19) 

N2(B)  ^  N2(A)-Hhv 

(20) 

N2(C)  N2(B)-thv 

(21) 

N2(A)  +  N2<r^  N2(v=0)+N2 

(22) 

N2(B)  +  N2  N2(v=0)-hN2 

(23) 

N2(a)  +N2  N2(B)-i-N2 

(24) 

N2(a)  -i-N  N2(B)-i-N 

(25) 

N2(C)  +N2  0  N2(a)-i-N2 

(26) 

N+N+N2  ^  N2(B)-I-N2 

(27) 

N-i-N-f-N  N2(B)-i-N 

(28) 

N-i-N-tN2  N2(A)-tN2 

(29) 

N+N+N  N2(A)-i-N 

(30) 

N2++  N  N2(v=0)-i-N+ 

(31) 

N-t  N  N2++e- 

(32) 

N+-I-  e”  N+hv 

(33) 

N2(a)  -i-N2(A)  N2(v=0)-i-N2+-i-e- 

(34) 

N2(a)  +N2(a)  ^  N2(v=0)-tN2++e- 

(35) 

N2(a)  +N2(v>24)  N2(v=0)+N2++e- 

(36) 

These  processes  are  considered  mainly  in  cold  discharge  models  and  the  data  are  quite  scanty  [11- 
17],  In  general,  for  those  processes  for  which  the  analytical  dependence  of  the  rate  coefficient  on  the 
temperature  is  unknown,  the  exothermic  rate  has  been  considered  constant  and  the  rate  of  the  reverse 
process  is  calculated  applying  the  detailed  balance  principle  [1-3]  which  assures  the  existence  of  the 
equilibrium  state. 
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This  reaction  scheme  has  been  improved  by  adding  electron-molecule  collision  processes,  in 
particular  vibrational  and  electronic  inelastic  (eq.  8)  and  superelastic  (eq.  9)  collisions.  The  rate 
coefficients  of  these  processes  are  calculated  from  the  eedf 

oo 

^vw  =  Jn{e)ovw(e)v(e)de 

0  (37) 

The  same  can  be  written  for  reactive  e-M  processes  such  as  dissociation  and  ionization.  The  cross 
sections  (o)  of  e-M  processes  can  be  found  in  [19-26].  It  should  be  noted  that  in  the  present  paper 
recombination  of  atomic  nitrogen  pumps  either  the  ground  state  through  the  ladder  climbing  model  or 
the  electronically  excited  states  through  reactions  27-30. 


RESULTS  AND  DISCUSSION 

The  model  previously  described  has  been  applied  to  the  F4  nozzle  [2]  (see  fig.  2). 


In  this  paper  we  investigate  the  role  of  electronically  excited  states  and  free  electron  kinetics  on  the 
F4  nozzle  flow.  At  this  purpose  we  have  compared  the  results  obtained  with  three  different  kinetic 
models; 

(a)  =  only  vibrational  and  vibrational-dissociation  kinetics  (proc.  10-rl5), 

(b)  =  (a)  +  electronically  excited  state  kinetics  (proc  10-;-36), 

(c)  =  (b)  -I-  free  electron  kinetics  (proc  8-r36). 

To  study  the  role  of  electron  kinetics  on  the  state-to-state  kinetics,  we  have  chosen  as  stagnation 
conditions  Tq  =  10000  K  and  Pg  =  1  atm.  In  previous  works  [5,6],  the  ionization  kinetics  have  been 
neglected,  supposing  constant  electron  molar  fraction.  In  reference  [7],  where  the  self-consistent  model 
includes  ionization/recombination  kinetics  for  the  conic  nozzle,  we  have  shown  that  electron  density  is 
practically  constant.  The  same  results  has  been  obtain  in  this  work.  On  the  contrary,  electronically 
excited  states  kinetics  [7]  seems  to  have  strong  effects  also  on  macroscopic  quantities,  as  can  be 
observed  in  fig.  3  where  the  mach  number  calculated  in  the  three  cases  a,  b,  c  has  been  compared. 
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x(m) 

FIGURE  3.  Mach  number  profile  calculated  in  different  kinetic  models;  a  proc  10-15,  b  proc  10-36,  c  proc  8-36. 


x(m)  x(m) 

FIGURE  4.  N2(A)  and  N2(a')  molar  fraction  profiles  calculated  for  different  kinetic  models;  a  proc  10-15,  b  proc  10-36, 
c  proc  8-36. 


Inclusion  of  processes  16-36  in  the  kinetics  has  dramatic  effects  not  only  on  the  electronically 
excited  states  (see  fig  4),  but  also  on  the  vibrational  distributions  (see  fig.  5)  mainly  for  vibrational 
quantum  numbers  v>10.  The  molar  fraction  of  N2(A)  is  influenced  by  e-M  collisions  only  in  the  region 
close  to  the  throat  (0<x<l  m).  When  e-M  collisions  are  neglected  (fig  4,  curve  b),  the  fraction  of  N2(A) 
is  larger  than  in  the  case  in  which  they  are  considered  (fig  4,  curve  c).  In  the  first  case  (b)  the 
recombination  of  atoms  in  the  state  N2(A)  produces  a  peak  around  0.2  m,  that  disappears  when  e-M 
collisions  are  considered  (case  c)  because  electronic  superelastic  collisions  dominates  over  all  the 
processes.  For  x>l  m  the  molar  fraction  is  not  affected  by  the  presence  of  e-M  collisions  because  the 
dominant  processes  are  the  quenching  of  N2(A)  by  atoms  and  molecules  and  atom  recombination.  The 
state  N2(a')  behaves  differently;  e-M  collisions  induce  higher  population  than  in  case  b,  due  to  inelastic 
collisions  in  the  region  close  to  the  throat  and  the  differences  remain  constant  for  x>l  m. 
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These  behaviours  are  reflected  on  the  vibrational  distributions  (see  fig  5).  Close  to  the  throat  (x=0.1 
m),  e-M  collisions  dominate  (in  particular  inelastic  e-v)  for  v<10  and  therefore  the  vibrational 
distribution  in  this  case  is  more  populated  than  in  the  other  cases,  on  the  contrary  the  distribution  tail  is 
determined  by  the  competition  betweenthe  atom  recombination  and  the  vibrational-electronic  coupling 
processes.  At  the  nozzle  exit  (x=4  m)  the  distribution  tails  are  populated  by  the  recombination,  while 
for  v<10  the  superelastic  vibrational  collisions  reduce  the  distribution  energy. 


FIGURE  5.  Vibrational  distributions  in  two  different  nozzle  position  calculated  in  different  kinetic  models; 
a  proc  10-15,  b  proc  10-36,  c  proc  8-36. 


FIGURE  6.  Vibrational  distributions  along  the  nozzle  calculated  for  model  c  (proc  8-36). 


To  complete  the  analysis  of  the  system,  in  figure  6  vibrational  and  electron  distributions  are  reported 
for  different  positions.  The  two  distributions  present  different  relaxations.  Vibrational  distribution 
cools  down  for  x<0. 1  m,  and  for  x>0. 1  m  the  recombination  becomes  very  efficient  and  the  tails  are 
populated.  Past  the  throat  all  the  distribution  are  non  Boltzmann.  The  eedf  has  an  opposite  trend.  For 
x<0.1  m  eedf  relaxes  very  slowly,  and  only  for  x>0.2  m  the  relaxation  speeds  up.  Moreover,  the 
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distributions  are  almost  maxwellian  till  x=0.5  m,  and  only  at  long  distance  the  eedf  presents  some 
plateaux,  a  characteristic  sign  of  the  presence  of  superelastic  electronic  collisions. 


CONCLUSIONS 

In  this  work  we  have  shown  that  the  relaxation  of  the  vibrational  and  electron  distributions  depends 
on  the  kinetic  model  considered;  electronically  excited  states  are  strongly  coupled  to  vibrational  kinetics 
and  electron  kinetics.  For  this  reason  self-consistent  models  become  the  only  tool  to  properly  describe 
plasma  kinetics.  The  model  can  be  easily  extended  to  investigate  plasma  jets,  considering  the  effects  the 
electric  fields  produces  on  the  kinetics.  The  results  reported  in  the  present  paper  are  qualitatively 
similar  to  those  presented  in  other  works  [6,7],  but  present  differences  due  to  the  nozzle  geometry.  The 
model  should  be  extended  to  air  mixtures  for  application  to  hypersonic  wind  tunnels  and  to  collisional- 
radiative  models  of  atomic  species  to  investigate  the  expansion  of  a  plasma  produced  by  laser  ablation 
[27]. 
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